Efficient atomic clocks operated with several atomic ensembles 

J. Borregaard 1 and A. S. S0rensen 1 

1 QUANTOP, Niels Bohr Institute, University of Copenhagen, 

Blegdamsvej 17, DK-2100 Copenhagen , Denmark 

(Dated: April 25, 2013) 

Atomic clocks are typically operated by locking a local oscillator (LO) to a single atomic ensemble. 
In this article we demonstrate a scheme where the LO is locked to several atomic ensembles instead 
of one. This results in an exponential improvement compared to the conventional method and 
provides a stability of the clock scaling as (aN)~ m > 2 with N being the number of atoms in each of 
the m ensembles and a is a constant depending on the protocol being used to lock the LO. We find 
an improvement for N > 20 for clocks operated with the conventional Ramsey protocol and N > 7 
for clocks operated with a more efficient, adaptive measurement strategy. 
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Atomic clocks have been subject to extensive investi- 
gations due to the importance of very precise time mea- 
surements in a broad range of areas in physics. The reso- 
lution of current atomic clocks is limited by the quantum 
noise of the atoms, which leads to the standard quantum 
limit where the resolution scales as 1/y/N with N being 
the number of atoms [TJ [2] . Various ways of improving 
the resolution have been suggested such as using entan- 
gled states with reduced atomic noise [SHI] to push the 
resolution to the Heisenberg limit where it scales as 1/N 
[THT2]. Another approach to increasing the stability is 
to use optical atomic clocks where the higher operating 
frequency leads to an improved stability [T3HT7] . Since 
an atomic clock is typically operated through Ramsey 
spectroscopy [18] the resolution can however also be en- 
hanced by increasing the Ramsey time T giving rise to 
an improvement scaling as 1/y/T [TM2T] . For clocks with 
trapped atoms, where there are no other limitations, T 
becomes limited only by the decoherence in the system. 
In practice this decoherence often originates from the fre- 
quency oscillations of the local oscillator (LO) used to 
drive the atomic clock transition [2D]. Hence instead of 
the advanced ideas discussed above the stability can also 
be increased by simply devising methods to increase the 
Ramsey period by stabilizing the LO. 

In this letter we suggest a scheme where the frequency 
of the LO is locked to the atomic transition using several 
ensembles of atoms. This procedure allows increasing the 
Ramsey period each time another ensemble is used. As a 
result we find that the stability of the clock can increase 
exponentially with the number of ensembles. Fig. [Ha) 
illustrates the general idea behind the scheme. The feed- 
back of the first ensemble locks the frequency of the LO 
thus reducing the noise to the atomic noise of the en- 
semble. Having reduced the noise in the LO, the second 
ensemble can be operated with a longer Ramsey time and 
through a second feedback loop the noise of the LO can be 
further reduced as shown in the simulation in Fig. [lib) 
(details of the simulation are given later). The proce- 
dure can be extended to any number of ensembles and 
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FIG. 1. (Colour online) (a): Illustration of the idea behind 
the scheme of locking the LO using several ensembles (a) . The 
feedback of the first ensemble stabilizes the LO such that the 
second ensemble can be operated with a longer Ramsey time. 
The feedback from the second ensemble then further stabilizes 
the LO. (b): Numerical simulation of the noise of the LO when 
locked to between 1 and 3 ensembles. The data was simulated 
as described in the text for TV = 20 and T\ — O.I/7. The first 
feedback lowers the noise of the LO and makes the spectrum 
white even though the unlocked LO was assumed to be subject 
to 1// noise. The second and third feedbacks further lowers 
the noise of LO by a constant factor. 



for uncorrelated atoms the stability of the LO will scale 
as ^(-fTiN)^" 1 ^ 2 where m is the number of ensembles 
(each with N atoms) , 7 is a parameter characterizing the 
frequency fluctuations of the unlocked LO, and Tj is the 
smallest Ramsey time in the setup i.e. that of the first 
ensemble. Hence the scheme can provide an exponen- 
tial improvement in the stability with the total number 
of atoms. In order for the clock to be stable we need 
7T1 <C 1 and hence the protocol requires a minimum 
number of atoms to improve the performance. With the 
conventional Ramsey protocol we find that the scheme 
works for a minimum ensemble size of N = 20 atoms. To 
further optimize the performance of the scheme we study 
an adaptive measurement protocol for estimating the LO 
frequency offset [22], which extends the applicability of 
the scheme down to ensembles with only 7 atoms. This 



makes the scheme relevant for atomic clocks based on 
trapped ions, which are typically constructed with only 
a few ions e.g. the ion clock in Ref. [19] was constructed 
with exactly N = 7. 

We will now describe the locking of the LO to the 
atomic transition using only two atomic ensembles op- 
erated with different Ramsey times. The idea described 
here can then easily be extended to a scheme with several 
ensembles. We model an ensemble of N atoms as a col- 
lection of spin-1/2 particles and define the total angular 
momentum J = J2i=i &i where Si is the spin-vector of 
the «'th atom. We define the angular momentum oper- 
ators J x ,J y and J z in the usual way and initially the 
atoms are pumped to have (J) along the z-direction, 
(Jx) = (Jy) = 0. In Ramsey spectroscopy the atoms 
are illuminated by a near-resonant 7r/2-pulse from the 
LO, followed by the Ramsey time T of free evolution, 
and finally another near-resonant 7r/2-pulse is applied to 
the atoms. The Heisenberg evolution of the operators 
J x ,Jy, and J z is J\ = J x , J 2 = sin(£0)J y - cos(6(j))J z , 
and J 3 = cos(S(f>)J y + sm(5(f))J z where 5(f) = 5wT is the 
acquired phase of the LO relative to the atoms. At the 
end of the Ramsey sequence J3 is measured and used 
to make an estimate S(/) e = — arcsin( J3/J) of 5(f> where 
J = N/2. The feedback loop steers the frequency of 
the LO towards the atomic transition by applying a fre- 
quency correction of Awj = —a5(f> e /Ti to the LO where 
a sets the strength of the feedback loop and T, is the 
Ramsey time of the i'th ensemble. 

The first ensemble is operated with free evolution time 
Ti and we assume that the second ensemble is operated 
with free evolution time T 2 = nT\ where n is an integer. 
Assuming that negligible time is spent on everything but 
the free evolution periods T\ and T 2 the frequency offset 
of the LO between time tk-i = (k— l)Tj and tk — kT\ is 

5u(t) = 5u (t) + Awi(t fc _i) + Aw 3 (t,„), (1) 

where 5ujq (t) is the frequency fluctuation of the unlocked 
LO, Awi(£fc_i) is the sum of the frequency corrections 
applied up to time tk—i from the first ensemble and 
Aoj 2 {tsn) is the sum of the frequency corrections applied 
up to time t sn from the second ensemble (s is found by 
rounding (k — l)/n down to the nearest integer). The 
iterative equations for Awi(ife_i) and Auj 2 (t sn ) are then 

Awi(tfc-i) = Awi(t*H2) - aS^ith-J/T! (2) 

Aoj 2 (t sn ) = Aw 2 (i (s _ 1) „) - a6<j> e2 {t sn )/T 2 , (3) 

where 8(f) ei (tk-i) and 5<fi e2 (t sn ) are the estimated phases 
from the first and second ensemble at times tk-i and t sn 
respectively. Using Eq. Q we can write the phase of 
the LO relative to the atoms of the second ensemble just 
before the measurement at time t sn as 

Sfa(t m ) = J 5u(t an -t')dt' 



where A(j)\_ 1 -, n = J Q 2 Aw2(t( s _i) n )di / is the accumulated 
phase due to the feedback of the second ensemble and 



8<t>{t sn )= / 8Q(t sn -t')dt' 



5uj Q {t sn -t') + Auj l {t { ^ l)n -t')dt' (5) 



is the accumulated phase due to the frequency oscilla- 
tions of the LO when locked by the feedback of the first 
ensemble. For now we assume that T 2 >• T\ such that the 
feedback of the first ensemble has stabilized the LO but 
later we will relax this assumption. From Eqs. ©-pi 
we then derive the difference equation 



h(t sn ) - S4>2(t( a 



84>(t sn ) - 0(t(»-i) n ) 
-a5(j) e2 {t(^_ 1)n ). 



(6) 



From this expression we see that the evolution of the 
second phase S4> 2 is essentially driven by the noise of the 
stabilized LO from the first step 5<j> but is stabilized by 
the second feedback loop described by aS(f> e2 . 

To solve Eq. ffib we need to characterize the width of 
the noise of the stabilized LO from the first stage, (Sep 2 ) = 
J^ 2 dtf^ 2 dt'(SCb(t)SCj{t')). From Eq. |J2) and (J5) we can 

derive a difference equation for 5(f)' {ti) = J Q 1 Sui(U— t')dt', 
which is the acquired phase of the LO relative to the first 
ensemble between time U-\ and ti (we can neglect the 
feedback of the second ensemble since T 2 ^> Ti): 



l'{ti) - ty'(ti_i) = SMU) - SMU-i) 
-aS^iU^). 



(7) 
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Here Sifiofti) — L 1 Su>o(ti — t')dt' is the phase of the un- 
locked LO and t( s _i)„ < £j < t sn . To solve this equation 
we follow Ref. |23J where the locking of the LO to a sin- 
gle ensemble is described. Considering for now a LO 
subject to white noise we write 5(f)' (t) ~ 5uj(t)Ti = x±(t) 
and 5(po(t) ~ 5u)o(t)Ti = t\Zi(t) where e 2 = 7T1 and 
(zi(t)zi(t')} — T\8{t — t'). Here 7 is a parameter char- 
acterizing the frequency fluctuations of the unlocked LO 
and 5{t — t') is a delta like function with width T\ such 
that Xi<5(0) = 1. We assume that N » 1 and that we 
can treat J X T-Jy,Jz as Gaussian variables. This enables 
us to derive a differential equation from Eq. (I7f of the 
form: 

ii(t) = €\z{t) — aarcsin(_ftT(l + e^z^t)) sin(xi(t)) 

+Ke 2 cos(x 1 (t))), (8) 

where we have introduced two additional noise fields 
z 2 ,z 3 fulfilling (zi(t)zj(t')) = 5i t jTiS(t 



A(f) 



(2) 

(s-l)n 



84>{t s 



(4) 



t') to describe 
the atomic noise. Si.j is the Kronecker delta, K — (J z )/J, 
e 2 = AJ y /(J z ), and e 3 = AJ Z /(J Z ). Note that we have 
here used 5(j> ei = — arcsin( J3/J). After expanding the 
arcsin function to third order it is straight forward to 



solve Eq. ^ and get an expression for (xi(t)xi(t')) . In- 
serting this in the expression for (S(f> 2 ) gives 

(<50 2 ) = 7 (7 2 7?(1 " K) 2 (l/6 - X + 4/9 X 2 ) 

+4/{ 1 T l ) + 4(l-x))T 2l (9) 

where \ = \ a (Jz) / J ■ For atoms starting out in a coher- 
ent spin state K = 1,AJ Z = and AJ y /(J z ) = 1/-/AT. 
In this case we find that the width of the noise of the sta- 
bilized LO from the first ensemble is (8(j) 2 ) = T 2 /NTi = 
7T2. Furthermore this noise will have a white spec- 
trum for both white and 1/ / noise in the unlocked LO 
(see Fig. [TJb) and Rcf. 23 ). The second ensemble thus 
sees an effective white noise spectrum of the LO with 
7 = 1/{T X N). 

We now return to Eq. (pi. In line with the above 
arguments we define 8<p 2 ~ Sur(t)T2 = £ 2 (t) and 8<p(t) ~ 
8Cj(t)T 2 = hzi{t) where e 2 = 7X2 and {zi{t)zi{t')) = 
T 2 5(t - t') (X 2 <5(0) = 1). The stability of the clock after 
running for a time r 3> T 2 is then 



much lower than the single ensemble result in Eq. (11) 
(with T 2 — >T\) and we have an improved stability. The 
arguments leading to Eq. ( 12 ) can be generalized in a 



a 7 (r) = ((fc(r) ; -i 
1 



2\l/2 



ujtTo 



dt / dt'(x 2 (t)x 2 (t')) 
Jo 



1/2 



, (10) 



where 8ui(t) is the mean frequency offset and u is the 
frequency of the atomic transition. Following similar ar- 
guments as before we can derive and solve a differential 
equation similar to Eq. ffity for x 2 (t) and obtain an ex- 
pression for (x 2 (t)x 2 (t')). Inserting this into Eq. (10) 
and taking the limit of r 3> T 2 results in 



straight forward way to show that if the LO is locked 
to m ensembles of N atoms the stability of the clock 
is ct 7 (t) = v / (/3i//32) (m - 1) 7/(w 2 r)(7V 7 T 1 )- m /2 ( sincc the 

noise of the LO is white after locking it to the first en- 
semble we use /?2 = 03 = • • • = fim)- By continuing the 
procedure we thus improve the stability exponentially! 

In our analytical calculations above we have assumed 
N ~^> 1. To investigate the performance for smaller N 
we have simulated an atomic clock locked to between 1 
and 4 atomic ensembles each with atom numbers from 
N = 20 to 100. From the simulations we can generalize 
to the case where the LO is locked to m ensembles. We 
simulate the full quantum evolution of the atomic state 
through the Ramsey sequences and subsequent measure- 
ments and implement the feedback on the LO similar to 
the description in Eq. 11]) and above. The assumption of 
T 2 ^> T\ can be relaxed by applying a phase correction 
in the measurement of the ensembles (see supplemental 
material) . The lowest number of atoms for which the pro- 
tocol works is then N m i n ~ 2/(7X1) since we need each 
feedback to lower the noise by at least a factor of 2 to 
allow twice the Ramsey period at the next level. Increas- 
ing Xi will thus both improve the stability of the clock 
and lower the number of atoms needed for the scheme 
to work. Increasing T± will, however, also increase the 
width of the noise of the unlocked LO, e.g. for white 



noise the width of the noise 



of the unlocked LO is 



<7 7 (t) 



1 
Wo 



1 



tNT 2 



(11) 



Eq. (11) describes how the stability improves with T 2 
and N. The longest T 2 we can allow is determined by 
how well the LO is stabilized by the first ensemble as 
contained in 7 and we parametrize it by T 2t7nax = ^2/7. 
In a similar fashion we assume that T\ max — fti/j for 
the first ensemble. With these parameterizations we can 
express the stability as 



7 



1 



(S(J)q) — Pi. For experiments or simulations running with 
a fixed Ramsey time there is always a finite probability 
that phase jumps large enough to spoil the measurement 
strategy occurs since Ramsey spectroscopy with projec- 
tive measurements is only effective for phases < tt/2. In 
our simulations we see these phase jumps as an abrupt 
break down as we increase the Ramsey time Xi. Simu- 
lating a clock with I = 10 6 results in a maximal Xi of 
fix ~ 0.1. This means that the scheme of locking the LO 
to several ensembles will only work for TV > 20. Ideally 
we should include correction strategies for the errors that 
limits Xi but for simplicity we ignore this. This means 
that our simulations have a weak (logarithmic) depen- 
dence on the number of steps that we simulate, which we 
do not expect to change our results significantly [2"2"| . 

We have simulated clocks with an unlocked LO subject 
to both white and 1// noise. The noise spectrum of the 
unlocked LO is defined as S(f)S(f + f) = (8u(f)8w(f')) 
where 8tu(f) is the Fourier transform of the frequency 
fluctuations 8ui(t) of the LO. For a free running LO with 
white noise we use S(f) = 7 while for 1/f noise we 
use S(f) = 7 2 //- The stability of the clocks is plot- 
ted against the ensemble size N in Fig. [2j Fig. [2] shows 
that the scheme works down to atom numbers of N = 20 
where we gain a factor of ~ 2 m ~ 1 in <t 2 (t) by locking 
If 7V7X1 3> 1 the stability obtained from Eq. (12) is the LO to m ensembles instead of one. This result holds 



cr 7 (r) = - 

LO 



-^V 2 /3i/3 2 



7/V& 



(12) 



Note that with the assumption of white noise in the un- 
locked LO we can pick j3i = fi 2 . As previously noted the 
noise of the LO will also be approximately white with 
7 ~ 1/iVXi after locking it to the first ensemble also for 
other types of noise e.g. 1/f noise. In that case it is 
desirable to have j3 2 ^ ji\ but we still expect /3i//3 2 to be 
of the order unity. Eq. ( J12] ) shows that by locking the 
LO to two ensembles of uncorrelated atoms the stabil- 
ity can be significantly improved from the performance 
obtainable by locking the LO to only a single ensemble. 
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FIG. 2. (Colour online) The stability of atomic clocks for 
a LO subject to (a) white noise and (b) l//-noise. • ■ A 
and T is the stability of a clock with the LO locked to 1,2,3 
and 4 ensembles containing N atoms each. The clocks were 
simulated with /3i = 0.1 and Tj = nTj_i where n = 0.1N for 
j > 1. The dashed lines are the analytical calculations. 



for both an unlocked LO subject to white noise and 1// 
noise. Furthermore the numerical results are seen to 
agree nicely with the analytical calculations. 

The conventional Ramsey protocol considered so far 
has a lower limit of N min ~ 20 in order for the scheme 
of locking the LO to several ensembles to improve the 
stability. This limit is due to the conventional protocols 
inability to effectively resolve the phase of the LO rel- 
ative to the atoms for phases > ir/2. In Ref. [22] we 
presented an adaptive protocol for estimating the phase, 
which effectively resolves phases < 7r. For a clock run- 
ning with I = 10 6 this protocol enables us to extend the 
Ramsey time to /3i ~ 0.3 if the phase fluctuations of the 
unlocked LO have a white spectrum and to /3i ~ 0.2 if 
the phase fluctuations have a 1// spectrum. The basic 
idea behind the adaptive protocol is to use a series of n 
weak measurements with intermediate feedback to rotate 
the atomic state to be almost in phase with the LO. The 
final estimate of the drifted phase is then the sum of the 
n rotations. Due to the rotations the protocol is able to 
distinguish whether the phase lies in the intervals [0; tt/2] 
or [tt/2; it] and thus avoids the ambiguity between these 
intervals in the conventional protocol. We assume that 
the weak measurements are performed according to the 
strategy developed and demonstrated in Ref. [9l QUI [24] 
where a light field dispersively interact with the atomic 
spin and is subsequently measured. With this adaptive 
measurement strategy we simulate clocks locked to be- 
tween 1 and 4 ensembles for atom numbers from N = 7 
to 34 with an unlocked LO subject to both white and 
1// noise. The stability of the clocks is plotted against 
the ensemble size N in Fig. [3] For the adaptive protocol 
we can apply the scheme of locking to several ensem- 
bles down to ensemble sizes of N = 7 where we gain a 
factor of ~ 2 m_1 in ct;?(t) by locking the LO to m ensem- 
bles instead of a single ensemble. This limit holds for an 
unlocked LO subject to both white and 1// noise even 
though Pi needs to be a bit smaller for 1// noise than 
for white noise. This is because the feedback of the first 
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FIG. 3. (Colour online) The stability of atomic clocks with 
adaptive measurements [22] for a LO subject to (a) white 
noise and (b) l//-noise.», ■ A and T is the stability of a clock 
with the LO locked to a 1,2,3 and 4 ensembles of N atoms 
each. The clocks were simulated with the adaptive protocol 
of Ref. 22 that allows for /3i = 0.3 (a) and /3i = 0.2 (b). For 
both noise types Tj — nTj-i with n = 0.3iV for j > 1. The 
dashed lines are the analytical calculations. 



ensemble effectively reduces the noise of the LO to ap- 
proximately the same level as for a LO with white noise 
operated with j3\ ~ 0.2. We can thus allow for f3 ~ 0.3 in 
the later ensembles since the noise of the stabilized oscil- 
lator is white. This means that the minimum ensemble 
size for which the scheme works is the same for 1// noise 
and white noise in the unlocked LO. 

In conclusion we have demonstrated a scheme for lock- 
ing the LO in an atomic clock to m ensembles of N atoms 
each. For this scheme the stability of the clock scales as 
v /7(7TiiV) _m / 2 where T\ is the Ramsey time of the first 
ensemble. Our scheme thus provide an exponential im- 
provement in the stability with the number of atoms. For 
the conventional Ramsey protocol our scheme is applica- 
ble down to ensemble sizes of TV = 20 atoms while it 
is applicable down to ensemble sizes of TV = 7 for the 
adaptive protocol considered in Ref. [22]. This make 
the scheme relevant for atomic clocks with trapped ions 
if the adaptive protocol can be implemented efficiently. 
An interesting continuation of this work is to investigate 
other adaptive measurements protocols, which might be 
easier to implement on trapped ions. Furthermore the 
performance of the protocol can most likely be enhanced 
by considering squeezed states but the treatment of this 
is beyond the scope of this article. 

During the final stage of our work a highly related 
preprint appeared, which treats essentially the same lock- 
ing scheme [25j . In that paper a different measurement 
procedure, different atomic states as well as dynamical 
decoupling is considered. Furthermore the upper limit 
to the Ramsey time is also viewed from a different per- 
spective. Taking this into account our results for the 
conventional Ramsey clock is more or less consistent with 
theirs but with the adaptive protocol we extend the ap- 
plicability of the scheme to ensembles of only 7 atoms. 
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SUPPLEMENTAL MATERIAL: EFFICIENT ATOMIC CLOCKS OPERATED WITH SEVERAL ATOMIC 

ENSEMBLES 

This supplemental material to our article "Efficient atomic clocks operated with several atomic ensembles" describes 
the details of our numerical simulations of atomic clocks locked to several atomic ensembles. We show how the 
assumption of T 2 ~^> T\ made in the article can be relaxed by applying a phase correction in the measurement of the 
second ensemble. We do this for the general case of locking the LO to m-ensembles. 

PHASE CORRECTIONS 

The Ramsey sequence and the subsequent estimate of the drifted phase of the LO relative to an ensemble of atoms 
is described in the article. Eq. (1) - (3) in the article describes the frequency offset of the LO (Sui(t)) between time 
th-i = (k—l)T and tk = kT when the LO is locked to two ensembles. We will now generalize this formalism to 
the case where the LO is locked to m ensembles. Assuming that the j'th ensemble is operated with Ramsey time 
Tj = n? T\ (n is an integer describing how many times the Ramsey time can be increased for each added ensemble) 
the frequency offset of the LO between time tk-i = {k— l)Ti and tk = kT\ is 

6u(t) = Suj (t) + Au x (t Sl ) + Auj 2 (t S2 n) + ■■■ + Aw TO (t, mn »-i), (13) 

where 5u>o(tk) is the frequency fluctuations of the unlocked LO and Auj m (t s . n j-i) is the sum of the frequency correc- 
tions applied up to time t s . n j-i from the j'th ensemble (sj is found by rounding (k— l)/n J ~ 1 down to the nearest 
integer). Note that the index SjH^ 1 should be read as Sj times v?" 1 describing the exponential increase in the 
Ramsey time each time another ensemble is used. The iterative equation for AWj(t s .„j-i) is 

Awj-(t 4jB i-i) = &Uj(t isi _ l)n i-i) - aS^(t Sin i-i)/Tj, (14) 

where S(j) ej (t s . n j-i) is the estimated phase from the j'th ensemble at time t s . n j-i and a sets the strength of the 
feedback loop (we assume equal strengths for all feedback loops) . a determines how long time the clocks needs to run 
before the LO is effectively locked by the feedbacks. In the article we assumed that T 2 ^> T\ such that the feedback of 
the first ensemble had effectively locked the LO before the measurement of the second ensemble. In the general setup 
of locking the LO to m ensembles this corresponds to assuming that n>l. We will now show how we can apply a 
phase correction in the measurement of the j'th ensemble such we can relax this assumption. The phase correction 
will compensate for the fact that the information from the last measurements on the first (j — 1) ensembles has not 
been fully exploited by the feedback loops before the measurement of the j'th ensemble. 

The phase of the LO relative to the j'th ensemble just before the measurement at time t sn j-i is 

(15) 
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(16) 

where 5<j> B +(s~x)ni- x IS the accumulated phase between time ti s . _i\ n j-i +s _ 1 and t( g -_iw-i +s due to the frequency 
fluctuations of the unlocked LO and the feedback corrections applied up to time t( s _iw-i. For simplicity we have 
replaced the time dependence by an index such that <5</5 e * i_ t , , ,w_, is the phase estimate from the j'th ensemble at 
time £ s /„i-i_|_( s ._i)„j-i. To fully exploit all information from the measurements on the first (j — 1) ensembles between 
time i( S3 _i) n j-i and t Sjn j-i, we choose a phase correction of ^ c °"°'ii : ' = </)™ 1 rrcct + (^correct j t 2 + . . . + ft'jj-i* where 

^correct,, ,; = £ (1 - a) ni ~ t -'S<t^ li . 1+(tf . 1)ni . 1 +c«^(l- o) Bi "'"'^? n .- 1+(ar .i )nJ - l (17) 

s=l s'=l 



For this choice of <J> corl ., C -i 3 the phase of the LO relative to ensemble j is 
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accumulated errors between the estimated phases and the actual phases for the (j — l)'th ensemble between times 
t(s-i)ni- 1 andt s . n .j-i (this is seen by considering Eq. (18) for j = 1,2, . ..). 3? J . n3 --i is thus the accumulated phase of 
the LO between time t(|._iw-i and t s . n j-i when it is fully stabilized by the feedbacks of the first (J — 1) ensembles. 
With the phase corrections Q^^-i' 1 we can therefore relax the assumption of n ^S> 1. As opposed to the feedback 
loop, which correct for e.g. frequency drifts by changing the frequency of the LO, the phase corrections directly 
correct the phase. This phase locking ensures a more rapid convergence, which is important when we want to apply 
the LO to the subsequent ensembles. Since the noise of the LO is white after stabilizing it to the first ensemble 
the subsequent frequency corrections from the other ensembles could be replaced with merely phase corrections of 
the LO, which would simplify the above procedure of applying a phase correction in the measurement of the j'th 
ensemble. We have however chosen to consider frequency corrections to keep a consistent treatment of the feedback 
in all stages. 



In our simulations we are simulating a clock with a LO locked to m ensembles running for a long but finite time. 
Similar to our description of the phase corrections $ corl j cx ; i ■? above there will be some remaining information from 
the last measurements, which have not been fully exploited by the feedback loops when our simulation stops. In 
our simulations we therefore do an additional phase correction $|° r ai cct to the LO after the final measurement. In 
principle the influence of the last few measurements could also have been reduced by running the simulation for a 
longer time but by doing the phase correction we reduce the required simulation time. With the phase correction the 
mean frequency offset of the LO (u>(t)) after running the clock for a total time of r = ZTi is 
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where 4> s = L 1 Su(t s — t')dt' is the phase of the LO relative to the atoms at time t s and ^^lf ct is the final phase 
correction of the LO. Using Eq. (13l-(14) and assuming that the j'th ensemble is operated with Ramsey time 
Tj = n? ' T\ we can write w(r) as: 
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(20) 



where 5<fi® is the accumulated phase between time i s _i and t s due to the frequency fluctuations of the unlocked LO 
and Srff 3 , i _ 1 is the estimated phase from the j'th ensemble at time t s , n j-i. We find that the ideal performance is 
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(21) 



With this phase correction the mean frequency offset is 
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Using Eq. ( 18 1 we can write: 
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(24) 



where T m is the Ramsey time of the m'th ensemble, $™ is the accumulated phase of the stabilized LO relative to the 
atoms in the m'th ensemble at time t sn m-i and <\>\ m is the estimate of that phase. We use Eq. (24 1 to determine the 
stability of the clock, which is given by cr 7 (r) = ((Suj(t) / 'uj) 2 ) 1 I 2 . 



